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We present a comprehensive study of non-equilibrium phenomena in the low temperature phase 
of the Edwards- Anderson Gaussian spin glass in 3 and 4 spatial dimensions. Many effects can be 
understood in terms of a time dependent coherence length, It, such that length scales smaller that 
It are equilibrated, whereas larger length scales are essentially frozen. The time and temperature 
dependence of It is found to be compatible with critical power-law dynamical scaling for small 
times/high temperatures, crossing over to an activated logarithmic growth for longer times/lower 
temperatures, in agreement with recent experimental results. The activated regime is governed by a 
'barrier exponent' ip which we estimate to be ip ~ 1.0 and ip ~ 2.3 in 3 and 4 dimensions, respectively. 
We observe for the first time the rejuvenation and memory effects in the four dimensional sample, 
which, we argue, is unrelated to 'temperature chaos'. Our discussion in terms of length scales allows 
us to address several experimentally relevant issues, such as super-aging versus sub-aging effects, 
the role of a finite cooling rate, or the so-called Kovacs effect. 

PACS numbers: 05.70.Ln, 75.10.Nr, 75.40.Mg 



But it's better to have little free time, then memo- 
ries don't intrude. Still, my God, what an amazing phe- 
nomenon these memories are. 
D. Shostakovich (Letter to I. Glikman). 

I. INTRODUCTION 

Spin glasses represent a model system for more 'com- 
plex' glassy materials. As sujch, it has attracted a large 
attention in the last decadesa. Despite the large number 
of theoretical, numerical and experimental papers pub- 
lished in the field, the original spin glasSpproblem is still 
far from being quantitatively understooda. 

Experimental facts that have to be explained are 
mostly of dynamical nature. The low-temperature dy- 
namics of various spin glass systems has been thoroughly 
investigated!™, and a number of characteristic features 
have emerged. Among these are the well-known aging 
behaviour of the thermoremanent magnetisation or the 
a.c. susceptibility in isothermal experimentsd, and the 
spectacular 'memory' and 'rejuvenation' effects observed 
in temperature-shifts protocols!!!]. 

It is fair to say that none of the existing theories pro- 
vides a complete quantitative ^description of spin glass 
dynamicsEI In a recent paperEJ, following earlier .ideas 
put forward in the context of the droplet modefcMJ, 
the notion of separation of length and time scales was 
argued to be crucial to account for the whole of exper- 
imental data. An important ingredient is the existence 
of a time and temperature dependent coherence length, 
^T(t w )y such that length scales smaller than £t are equi- 
librated, whereas length scales larger than £x are frozen. 
As a resultO~ll3, the aging dynamics after a certain age 



t w is ascribed to the motion of objects of size ~ £r(tw)- 
Note that we do not need to specify the topological na- 
ture o£-these objects, postulated to be compact droplets 
in Rcpli-but that could be also fractal, sponge- like struc- 
turesEjj~t2l. A simple realization of the above scenario 
was recently worked out in the context of the 2D XY 
modeJlj, suggesting that systems with quasi-long range 
orderEj have a dynamics very similar to what is observed 
in spin glass experiments. 
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FIG. 1. Growth law of the coherence length, Eq. (|l|), for 
T/T c = 1.0, 0.9, •■• 0.5 (from bottom to too). The numer- 
ical values are extracted from experimentstJ. The separa- 
tion of length scales within a given time window is very clear 
from these curves. Note also the crossover from the critical 
(power law) regime to the activated regime, where the coher- 
ence length is stuck to a value of a few £(T). 
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In this paper, following some previous works00 _ @, 
we identify a coherence length that grows with the ag&-af 
the system in the Edwards- Anderson spin glass modelEa. 
Wc .numerically test the growth law of £t proposed in 
Rcf.y, according to which: 
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where z c is the dynamical critical exponent, ip the so- 
called barrier exponent that describes the growth of the 
energy barriers with length scales, and T(T) a temper- 
ature dependent free-energy scale that vanishes at the 
critical temperature T c . This growth law, illustrated in 
Fig. 0, is motivated both-by theoretical considerationscJ 
and experimental results L1L and was used to analyze fur- 
ther experimental dataoEO. 

By construction, this growth law reduces to usual criti- 
cal scaling when the barriers on scale £t are much smaller 
than T. Assuming that the barrier scale behaves as 
T(T) = XJ1~T/T C )^ (where v is the correlation length 
exponential, and To an energy scale of order T c ), the 
crossover between critical scaling and activated scaling 
occurs for a dynamical crossover length £(T) that di- 
verges at T c as (1 - T/T c )- V . Remark that Eq. (Q) is ob- 
viously not the only possibility to describe this crossover. 
However, this multiplicative form was found to represent 
quite accurately, the .dynamics of the directed polymer or 
the Sinai modelr 8 r 9 l, where a similar crosscttier between 
diffusive and activated dynamics takes placeEH. 

As noted in RefO, the growth law Eq. (|l|) is difficult 
to distinguish, over a restricted range of length scales, 
from a pure power law t w ~ tq£ 7 t with a temperature 
dependent exponent z = z(T) > z_.|-Tlie latter was 
previousl*i-|reported both numericallyoE3~t3 and exper- 
imentalljfii. However, we believe that Eq. (|l|) should 
be prefered. One reason is that more elaborated experi- 
mental protocols, such as temperature-shift experiments, 
reveal non-activated effects, as recalled below. This non- 
activated behaviour is captured by Eq. (|l|), both through 
the temperature dependent barrier term T(T) and the 
strong renormalisatiDEL.aLthe microscopic time scale by 
critical fluctuation^E'El 

The present work is a quantitative investigation of 
the low-temperature, non-equilibrium, dynamics of the 
Edwards- Anderson spin glass modeled in finite dimen- 
sions, d = 3 and d = 4.ptt can be viewed as the numer- 
ical counterpart of RefO. Here, we take advantage of 
the fact that simulations, unlike experiments^, directly 
give access to It, to confirm some of the results obtained 
in Ref.til using indirect evidence. To do so, we perform 
an extensive series of numerical experiments, including 
simple aging, temperature-shift and temperature-cycling 
protocols. We observe for the first time in this system the 
'rejuvenation and memory' and 'Kovacs' effects, which 
are interpreted using the coherence length ^r(iiu). In 
turn, this allows us to shed new light on several ques- 
tions such as sub-aging effects, the issue of temperature 



chaos and the existence of an overlap length, and the very 
nature of the spin-glass phase. We emphasize also that 
although simulations and experiments are performed on 
very different time windows, the length scales probed dy- 
namically are actually not very different, see Fig. [j]. 

The paper is organized as follows. Section y intro- 
duces the model and gives technical details on the simu- 
lation. Section III focuses on simple (isothermal) aging. 
The growth law of the coherence length is studied in sec- 



IV. 'Small' temperature-shift experiments are per- 
while 'larger' shifts and cycles are 
Physical implications of our results 



tion 

formed in scctionl 
studied in section [ 
are discussed in section VII , and section VIII summarizes 
and concludes the paper. 



II. MODEL AND TECHNICAL DETAILS 

We study the Edwards^Anderson spin glass model de- 
fined by the HamiltonianCJ 



Hj[s] 



(2) 



where s — {si}i=i,... ,n are N = L d Ising spins located 
on a 3d or Ad (hyper)cubic lattice of linear size L, and Jij 
are random variables taken from a Gaussian distribution 
of mean and variance 1 . The sum is ovepjjiearest neigh- 
bors. The spin glass transition is believecO to take place 
at T c (d = 3) = 0.95 and T c (d = 4) = 1.8. In all this 
paper, the temperature is given in units of the critical 
temperature, T/T c (d) — ► T. 

To study the aging dynamics, we use a rather large sys- 
tem linear size L, 30 < L < 40 in d = 3, and 15 < L < 26 
in d — 4. On the time scale of the simulation, the system 
never equilibrates on a length scale larger than, say, ~ 8 
lattice spacings, and we thus always work in the regime 
(-r{t) L. The dynamics associated to the Hamiltonian 
(||) is a standard Monte Carlo algorithm, where the spins 
are randomly sequentially updated. One Monte Carlo 
step represents N attempts to update a spin. 

The behavior of the system is analyzed through the 
measurements of various physical quantities. 



We compute the energy density, defined by 



1 



e=-(Hj[s}), 



(3) 



stands for an average over initial con- 



where ( 

ditions and • • • over the disorder. 

• We measure the two-time autocorrelation function 
defined by 



1 



JY 



+ t,t w ) = ( Si (t w +t) Si (t w )). (4) 

i=l 
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We will also consider an ax. susceptibility like 
quantity x(w,tw), defined ascj: 



1 — C(t w + — , t w ) 



T 



(5) 



As in previous studies, we extract a coherence 
length by studying the dynamical 4-point correla- 
tion function, defined as 

1 N 

(6) 



i=l 



where (a, b) are two copies of the system starting 
from different initial conditions and evolving with 
independent thermal histories. This four-point cor- 
relation function can be interpreted as the proba- 
bility that two spins separated by a distance r have 
the same relative orientation in two independent 
systems after time t w , as is measured by a two- 
point function in a pure ferromagneto. 

Our data are typically averaged over 15 (autocorrela- 
tion functions, energy density) to 50 (4-point correlation 
functions) realizations of the disorder. Our data are thus 
reported without errorbars, which are typically extremely 
small. Technically, the most difficult part of our work is 
then to perform a meaningful analysis of the data in or- 
der to extract quantitative values of the various physical 
parameters. 



through careful noise measurements^. It is represented 
as a function of the time difference t in Fig. 2] for the 
4 dimensional sample. Similar curves are obtained at 
all temperatures, in d = 3, 4. We get the well-known 
'two-step' decay of the correlation function with a first, 
stationary, part followed by a second, non-equilibrium, 
aging part. The existence of these different regimes is 
easily understood qualitatively, but a more quantitative 
description of both time sectors is not completely settled 
yet. 




FIG. 2. Two-time autocorrelation functions C(t w +t,t w ) 
in simple aging experiments in d — 4 and T = 0.7. We show 
20 different waiting times, which are logarithmically spaced 
in the range t w £ [2, 57797] and increase from left to right. 



III. ISOTHERMAL AGING: BASIC FACTS 

In this section, we consider isothermal aging proto- 
cols. The system is quenched at initial time t w — from 
an infinite temperature to a low temperature T < 1.0 
where it slowly evolves towards its equilibrium state. 
Although the phenomenology is very-well knownB'lIiJ, 
and h as .. alre a dy been thoroughly investigated in simu- 
lationsoE3nJH, some important points are still poorly 
understood. We discuss all these aspects in some details 
in this section. We evaluate, in particular, the impli- 
cations of Eq. (1) for the theoretical description of the 
data. 



A. The spin-spin correlation function: general 
considerations 

As is now well-documentedlll, the slow evolution of 
glassy materials following a quench is best analysed 
through the measurement of a two-time quantity, typ- 
ically susceptibility or correlation functions. Here, we 
measure in the process of isothermal aging the two-time 
spin-spin correlation function C(t w + t,t w ) of the sys- 
tem. This quantity can also be accessed experimentally 



Ekjom a theoretical point of view, both mean-field mod- 
elstia and the multi-layer trap modeKl predict that the 
short-time and long-time contributions are additive, 

C(t w + t, t w j — C C q(t) + C a gi ng (t, t w ), (7) 

whereas-aging at a critical point leads to multiplicative 
scalingsEH, 

C(t w + t, t w ) — C C q(i)C a gi n g(t, t w ). (8) 

as used both in pimulationsEi and in early analysis of ex- 
perimental datauH. The equilibrium part can be fitted, 
both experimentally and numerically, by a power-law 

C eq (t) ~ Ar< T \ (9) 

with a temperature dependent exponent x(T), which 
takes rather small values. These two forms (additive and 
multiplicative) are actually not very different for short 
times, since C ag i ng (f, t w ) is approximately constant for 
t <C t w , in the regime where C eq (t) varies most. How- 
ever, one should stress that the extrapolation of the mul- 
tiplicative scaling behavior (Q) associated to (^|) to large 
times implies a zero Edwards- Anderson parameter, de- 
fined dynamically as: 

<7ea = hm lim C(t w + t, t w ). (10) 
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Indeed, no clear plateau appears in the curves of Fig. 
On the other hand, the additive scaling suggests a non- 
zeropvalue of qea, and accounts well for the experimental 
datsHH 

Various scaling forms have been predicted for the ag- 
ing contribution. In mean-field models, one expects an 
'ultramctric' bchavioi£3 



with the waiting time as 



Caging (t: ^w) ^ ^ 



hj(t w + 1) 
hi(t w ) 



(11) 



where the infi.nite.sum over the index H' refers to various 



unctions Ci and hi 
An explicit exam- 



'time sectorsOO, and the various 
have yet unknown functional formal . 
pie of such a scaling has recently been given in Rcf.EJ, 
in the context of the trap model, where the infinite sum 
boils down to 



Caging (^s tw ) C 



In t, n 



(12) 



Note that it is \nt and not ln(t w + t) that appeats4ii 
this equation, which ensures dynamic ultrametricityr 9 lr°l. 
This scaling is not observed experimentally, except per- 
haps for t -C t w (see RefJ23). The scaling ( |l2| ) is similar 
to, but different fromjj-th^scaling form suggested by the 
droplet picture, wheraLj'EJ, 



Caging (t^w) — C 



i T [tw) 



(13) 



m.thl T {t w ) ~ (ln^) 1 /^. The 'droplet' scaling variable 
ln(t + t w )/hxt w suggests super-aging, i.e. an effective 
relaxation time growing faster than the age of the system 
t w , which is not borne out by experimental data showing 
instead a tendency towards sub-aging. We come back to 
this point below. 

In the absence of any compelling theoretical descrip- 
tion, both experiments and simulations have been phe- 
nomenologically fitted with some scaling functions of the 
type 



Caging \pi t"w) — C 



h(t + t v 
h(t w ) 




(14) 

unctional 
" issue of 
ely used 
h{u) 



where— .the function h{u) is given var, 
formal3, related to the often debatedull 
sub-aging versus super-aging behavior. A, 
form for h(u), which we adopt here, ii 
exp [u 1_M /(l — fi)] , where the exponent \i allows one to 
interpolate between super-aging (/i > 1) and sub- aging 
(fi < 1), via simple aging (jj, = 1, for which h(u) = u). 
The effective relaxation time is indeed given by i ro i ~ t^. 
Note that if one takes h(u) ~ £t(u) with £t given by 
Eq. (|IJ) , then, as in the droplet model, super-aging would 
also be observed for long waiting times, since the effec- 
tive relaxation time, defined as h(t w ) / h' (t w ) , now grows 



^rel (tw ) ^ t h 



Z c + 



T(T)V4^ 
T 



(15) 



Note finally that the existence of a growing coherence 
length iritw) in spin-glasses does not necessarily implies 
that the correlation function can be expressed in terms 
of this length scale only. Indeed, since time scales are 
broadly distributed, processes corresponding to different 
lengths scales are expected to-onix together. This may 
also happen in simpler modelsEi. 



B. Numerical results 

We now discuss our numerical results. We first show 
in the inset of Fig. |3| (top) that the simple scaling C(t + 
tw,tw) ~ C(t/t w ) obviously fails in describing the data. 
Neither short nor long time scales are correctly described 
by this form. 
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FIG. 3. Rescaled autocorrelation functions according to 
the scaling form (§) in d = 3 (top, T = 0.6, x = 0.0128) 
and d = 4 (bottom). The inset of the top figure is the simpler 
scaling form C(t + t w ,t w ) = C(t/t w ). 

We then show in Fig. || that when the short-time dy- 
namics is taken into account through the multiplicative 
scaling form (Q), the collapse looks almost perfect in 
d = 3 (Fig. ||, top), whereas a small super- aging trend 
subsists in d = 4 (Fig. ||, bottom). Indeed, in this t/t w 
representation, older curves are still above the younger 
ones, suggesting that rescaling the time by t w is not suffi- 
cient to superimpose all the curves. Hence, the introduc- 
tion of another fitting parameter is required to describe 
the data, namely the exponent /i defined above. 
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FIG. 4. Top: Rescaled autocorrelation functions according 
to the scaling forms (^|) and (|l4|), thus allowing for fi 1. 
Bottom: Rescaled autocorrelation functions according to the 
scaling forms (Q) and (|l4|), choosing parameters so that data 
extrapolate to a non-zero Edwards-Anderson parameter, as 
shown by the horizontal dashed line. 

We find that with (x, fi) as free parameters and the 
multiplicative form (0), the data can be nicely collapsed 



for the whole temperature range studied, 0.4 < T < 1.0, 
in both dimensions d = 3 and d = 4. An example of 
such a rescaling is given in Fig. [| (top). In d = 3, our 
data are consistent with ji = 1 in the whole temperature 
range, and we find an exponent x(T) in close agreement 
with values reported in Refc3. In d p.4, our finding for 
x(T) also follow the reported valuescJ. In addition, as 
suggested by Fig. ||, we find that the exponent \i has to 
be larger than 1 for T < 0.6, while (J, = 1 is compatible 
with the data for 0.8 < T < 1.0. This observation was. 
never reported, although a re-analysis of published dataEj 
confirms this trended. In both dimensions, we find that 
the scaling function C(x) behaves as C{x) ~ const when 
x <C 1, and as C{x) ~ x~ x for x 3> 1, as in Ref.E3. 

Close to T c (i.e. T = 1.0), we find /i = 1. This is 
physically expected, since standard non-equilibrium crit- 
ical dynamics gives indeed the scaling (|8|), with h{t) ~ 
£x(t) ~ t 1 /*". In that case, the cohm-fince length is the 
usual dynamic correlation lengthEZrCJO. The fact that 
the scaling function h(u) changes when the temperature 
is lowered in d = 4 suggests that the dynamics leaves the 
critical regime. It is thus a priori surprising that the mul- 
tiplicative power-law scaling still holds for low tempera- 
tures. The situation appears different in d = 3, where 
the dynamics does not show any clear change when the 
temperature is lowered below T c . 

One can therefore try to rescale the 4d data accord- 
ing to an additive scaling, which allows for a non-zero 
Edwards- Anderson parameter, Eq. (]?]). In this case, we 
have three free parameters, (x, /i, A), where A is the am- 
plitude of the stationary part A. This is unfortunately 
too much since in this case A (and tkus Qea) is very 
poorly constrained. As noted in Ref.cHl, the values of 
the parameters A and x(T) are in fact strongly anti- 
correlated, but the data is insufficient to pin down their 
individual values, and hence to conclude on the value of 
qea- However, noting that C ag i n g(^) > 0, Vx allows us 
to give a possible range for (?ea- For instance, for d = 4 
and T = 0.6, we find that <7ea € [0.57, 0.68] leads to a 
reasonably good rescaling of the data. These values are 
in agreement with prnvious estimations of the Edwards- 
Anderson parameteifij. 

Interestingly, the additive procedure has little impact 
on the value of the aging exponent \i. For d = 3, we still 
find that fj, = 1 allows for a good rescaling of the data, 
whereas in d — 4 the super-aging tendency is slightly 
reinforced by this rescaling. An example of this is shown 
m Fig. | (bottom). 



C. Conclusion 

From the above analysis of our numerical results on 
isothermal aging, we conclude on the following. 

(i) Although data are compatible with the standard 
scenario where q^A > for T < T Cl we cannot rule out, 
from our numerical study of spin-spin correlation func- 
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tion, the fact that asymptotically <7ea = both in d = 3 
and in d = 4. Longer simulations in d = 4l3 however 
seem to favor the additive scaling over the multiplica- 
tive scaling, and therefore a non-zero Edwards- Anderson 
parameter. 

(ii) The exponent x that describes the short-time decay 
of the correlation function also describes the behavior of 
the equilibrium a.c. susceptibility. The measurement of 
the latter allows then to obtain the value of x indepen- 
dently. It was then shown experimentally tiat when this 
is done, the additive scaling (0) works welloH. 

(hi) Whatever the chosen rescaling for the short-time 
dynamics, we find a systematic super-aging behavior in 
d = 4. This is consistent with the identification of 
the scaling function h{u) with a coherence length £t(u), 
growing like Eq. (pi) Indeed, the resulting relaxation 
time given by Eq. ( fl5[ ) can be written approximately as 
tt with 



1 



1 

Zc, 



(16) 



On the other hand, no super-aging is found in d — 3, al- 
though, as shown below, Eq. (|I|) also seems to hold. How- 
ever the distinction between Eq. (|l|) and a pure power- 
law with a temperature dependent exponent will be much 
harder to establish in d = 3. 

(iv) Our data show no tendency towards sub-aging, 
in contrast to what is consistently found in all experi- 
mentso. We shall see below that the introduction of a 
finite cooling rate (instead of the direct T = oo — > T 
quench considered in this section) in fact results in an 
effective sub-aging. 



IV. GROWTH OF A COHERENCE LENGTH 

We now turn to a more geometric characterisation of 
aging, and try to associate the stationary part of the cor- 
relation function to equilibrated small scale dynamics, 
and the aging part of the correlation function to out- 
of-equilibrium, large scale dynamics. The time depen- 
dent crossover scale is a coherence length, that .ssipuld 
be the domain size in a coarsening ferromagnetic, 
the dynamic correlation length at the critical pointo. 
In the .-ease of spin-glasses, a more subtle definition is 
neededcil'E3. As for the autocorrelation function, we dis- 
cuss in detail the physics involved in the scaling form 
of this function. Moreover, we extend previous works in 
d = 4 to a larger temperature range. This is necessary in 
order to extract the parameters involved in Eq. ([!]) . The 
latter analysis is also new in d = 3, and allows a direct 
comparison with experiments. 



A. The four-point correlation function 



1. Definition 

As proposed by several authors0il, the coherence 
length can be measured in a simple isothermal aging pro- 
tocol from the spatial structure of the 4-point correlation 
function C4 (r, t) , defined in Eq. (||) . This 4-point func- 
tion is the analog of thje structure factor in a usual do- 
main growth problems, adapted to the case of the dis- 
ordered system under study, where any growing pattern 
is random and can only be identified by comparing two 
independent real replicas of the same system, prepared 
at t = in a different random state. It measures the sim- 
ilarity of the relative spin orientations in the two systems 
at a distance r after time t. Typical results are presented 
in Fig. m The spatial decay of this correlation function 
becomes slower when the time increases, clearly indicat- 
ing the growth of a length $eale is, the system. This was 
already noted several timestaOcj. 
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FIG. 5. 4-point correlation function for T — 0.7 in d — 4 
measured in an isothermal aging experiment. Times are log- 
arithmically spaced in the range t m G [2, 57797] and increase 
from left to right. The horizontal dotted line corresponds to 
9eAi using qea = 0.4 from static studies. 



2. Functional form of C4 

The correct identification of the coherence length 
(■T{tw) is however not completely straightforward. In- 
deed, the naive definition 



C 4 (r = £ T (t w ),t w ) = c, 



(17) 



where c is an arbitrary constant, say c = 0.1, leads to 
inconsistent results, because the decay of £L is not purely 
(or even possibly stretched) exponentialeS. This fact is 
very clear in d = 4, where the definition ( |l7| ) with c = 
0.1 leads to a coherence length which is such that for 
some t w , IrAtw) > ^r 2 (*w) when T\ < T2, i.e. a faster 
growth at lower temperatures. This result is physically 
unacceptable. 
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From a more careful analysis of the dataea, one finds 
that C4 receives two contributions: a 'quasi-equilibrated' 
decay for r < ^t(^), followed by a 'non-equilibrium' 
decay at large distances. This reflects the fact that at 
time t w , the system has equilibrated up to a length scale 
^r(itu), with a non trivial equilibrium correlation func- 
tion. This is the direct analog of the 'two-step' behavior 
observed in the autocorrelation function. 

As suggested by previous studiesE§E3, a possible func- 
tional form is 



C 4 (r, t w ) 



1 



C 4 



e T (t w ) 



(18) 



with a temperature dependent a(T), and Ca{x) a scaling 
function. It is difficult to confirm or dismiss this result, 
since the numerical correlation functions typically decay 
over 3—5 lattice spacings only, and other functional forms 
are possible (see below). Note that very few equilibrium 
data is available for this correlation functioned, which 
would be a very interesting information to compare with 

As noted in Refill, Eq. (|18|) suggests that C 4 (r -> 
00, t w — 00) tends to zero. This must be contrasted with 
the prediction of the droplet picture or any other theory 
in which the-pverlap distribution is a trivial (5-function at 
Qea, whereto 



lim lim C^r, t w ) 

r — >oo t w — >oo 



2 

<?EA 



0(r~ 



(19) 



where (?ea is the Edwards- Anderson parameter and 6 the 
energy exponent, estimated to be ~ 0.2 in d = 3li3 and ~ 
0.7 in d = 4EEI, i.e., smaller than the values of a reported 
in Table | below. 

Although Fig. H suggests that 64(7-, t w ) is rapidly much 
smaller than <z E Aj and compatible with C^(r — > 00, 00) = 
0, we find that it is still possible to rescale the data 
according to scaling forms which imply a non-vanishing 
large distance limit. In order to illustrate this point, we 
tried the following ansatz (not motivated by any theoret- 
ical argument): 



C 4 (r, t w ) = (<jf A + a exp(-br)) C\ 



1 I 



\lT(t w ) 



(20) 



which allows a rescaling of the data as good as Eq. (18). 
In this case, the stationary part of C4 indeed tends to- 
wards <z!aj but much faster than r~ e . More generally, the 
inequality a(T) > 6 makes Eq. ( |l9| ) rather unplausible. 



3. Discussion 

A word of caution is however needed here. Although a 
non-zero value of the Edwards- Anderson parameter q^A 
is expected in the spin-glass phase, dynamical evidence 
for this is still rather weak, at least in d = 3. For ex- 
ample, as discussed in the previous section, the dynamic 



spin-spin correlation function cannot rule out a non-zero 
value of qEA- It is well-known that on the time scale 
(resp. system size) of dynamic (static) simulations, the 
apparent value of qea constantly shifts towards with 
increasing times or sizescil. In that sense, the evidence 
that C^ir, t w ) tends to zero at large distances could be 
compatible with a very small value of ^ea- The evidence 
against the simplest droplet picture is however stronger 
in d = 4, since the numerical evidence for q^A > is 
more compelling in this case. 

An alternative interpretation in three dimensions is 
that 9ea = at all temperatures, which means that ei- 
ther there is no true spin glass transition, or else that 
the nature of the spin glass phase is different from what 
has insofar been theoretically expected. This issue might 
also be related to the existence of large excitations of 
finite energy recently found in Ref.Ej. For instance, 
one could be in a Kosterlitz-Thouless (kt) like situa- 
tion where <7ea = 0, but the 4-point correlation function 
changes from exponential for T > T c to power-law for 
T < T c with, possibly, a temperature dependent expo- 
nent a(T). Along this (speculative) line of thought, it 
has been pointed out recently that the dynamics of a 
critical phase (such as the kt phaaa) shares many simi- 
larities with spin glass dynamicalacJ. 



B. Numerical results for the coherence length 

We thus adopt a phenomenological definition of It fa) 
as the time-dependent length which leads, using Eq. fllq), 
to the best numerical collapse of C^r, t w ) measured at 
different times. 

It is important to note that the numerical value of the 
exponent a(T) used in this scaling plot has a significant 
influence on the resulting growth law for iritw)- Since 
the spatial support of C4 is very small, it is impossible 
to determine this exponent numerically with great accu- 
racy. The conclusion is that even using the above scaling 
procedure, there is still some degree of arbitrariness in 
the definition of the coherence length It- 

We first report in Table | the values of a(T) found in 
our simulation for the d = 3 and d — 4 cases. Our values 
for d — 3 are quite close to the ones found in Rcf.El 
For d = 4, only the_-sjalue of a(T — 1.0) = 1.63 has 
been reported in Ref.tij, in excellent agreement with our 
determination. 

An example of data collapse for the 4-point correlation 
function can be seen in Fig. ^. As in Rcfj23, we find that 
the cut-off function C4 is compatible with a 'stretched' 
exponential form, C±(x) ~ exp(— x^) 1 but with an expo- 
nent (3 > 1. 



T 


a(T) in d = 3 


a{T) in d = 4 


1.0 


0.6 


1.63 


0.9 


0.5 


1.35 


0.8 


0.5 


1.25 



7 



0.7 


0.5 


1.0 


0.6 


0.5 


0.9 


0.5 


0.45 


0.9 



TABLE I. Values of the exponent a(T) in d = 3 and d — 4. 
Note that a(T) is nearly constant in d = 3 but significantly 
changes with temperature for d = 4. This will turn out to be 
important in the following. 



The growth of £r(t w ) for different temperatures, 0.5 < 
T < 1.0, is reported in d = 3 and = 4 in Fig. 0. 
From Eq. ([!]) , we expect the coherence length to grow as 
a power law at short times. This is the critical regime 
characterized by iritw) <C £,(T) such that the growth 

law is £T(t w ) ~ tw* c - For larger times, the exponential 
activated term in Eq. ([l]) should slow down the growth 
of ^r(i) m a temperature dependent manner. The locus 
of the crossover itself must be temperature dependent. 
All this is observed numerically, see Fig. ^. This is also 
qualitatively consistent with the law extracted from ex- 
perimentstlftJ and reported in Fig. [0. 




FIG. 6. Rescaled 4-point correlation function for T = 0.7 
in d = 4. The full line is a fit of the scaling function C4(x) to 
a 'stretched' exponential form, with /3 — 1.7. 

We are now in position to compare the growth law 
obtained numerically to Eq. (|l]) . The critical exponent v 
is taken from previous numericaLwork. We take v = 1.65 
in d = 3 and v — 0.8 in d = 4c3. The exponent z c and 
the microscopic time To are fixed by the data at T = T c . 
We find z c ~ 7.0 and t ~ 2.0 in d = 3, z c ~ 5.9 and 
To ~ 2.2 in <i = 4. The values for the dynamic exponents 
are compatible with previous determinationsaEl 

We are thus left with tp and To as free parameters. We 
find that Eq. ([!]) accounts very well for the data in d = 4 
with ip ~ 2.3 and To ~ 0.6, see Fig. fj]. In 3 dimensions, 
we were not able to use Eq. (0) with a fixed tq. Instead, 
the fits reported in Fig. ^give ip ~ 1.0 and T ~ 5.5, but 
were obtained by letting To to be temperature dependent, 
with a non-monotonic temperature behaviour, for which 
there is of course no physical explanation. Simpler power 
law fits with a temperature dependent z(T) — z c T c /T 
but a constant to = 1 give equally good results, with 
much less free parameters. This might indicate: 

• cither that, as discussed above, the whole Ising 
spin-glass phase in three dimensions is Kosterlitz- 
Thouless like, where the dynamics is indeed de- 
scribed by power-laws for all temperatures. This 
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would be compatible with the fact that no cross- 
over beyond the critical regime is detected in the 
evolution of C(t,t w ) and C±{r, t w ) in d = 3; 

or that Eq. ([!]) is too simple to quantitatively re- 
produce the detailed crossover between critical and 
activated dynamics in d = 3__tt must be noted 
that, as emphasized in Refs.E^lH, the simulations 
are performed right in the regime where the influ- 
ence of the critical point is still strong; 



• or else that the procedure to extract £r{t v 
Ci(r,t w ) is somewhat biased. 
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FIG. 7. Growth laws of the coherence length in 3 dimen- 
sions (top) and 4 dimensions (bottom). The points are the 
data, the full lines are fits to Eq. (|IJ). 

We should nevertheless add the following remarks 
about the three dimensional case. 

(i) As will be clear below, the simple power-law growth 
of It with z(T) = z c T c /T cannot explain the small tem- 
perature shift effects, that suggest - both numerically 
and experimentally - deviations from a pure activated 
growth. 



(ii) One can also extract from the data the local slope 
of logt(^T) as a function of log^T, which should be in- 
dependent of £t for a pure power law. One finds instead 
systematic deviations, such that the effective exponent z 
indeed increases with £t, as predicted by Eq. (Q). More- 
over, the amplitude of these deviations vanish as T in- 
creases towards T c , in a way very much compatible with 
Eq- ©. 

(iii) Finally, it has been suggested in the past that a 
pure power-law behaviour for iritw) is associated with 
'Replica Symmetry Breaking', which predicts that the 
whole low temperature phase is 'critical'. We disagree 
with this point of view: the growth of £r(t w ) could 
be asymptotically logarithmic, as in the droplet picture, 
even if the equilibrium phase is not unique. This seems 
to be two totally separate issues as long as one does not 
associate It with the size of compact droplets. 



V. PROBING THE BARRIERS: SMALL 
TEMPERATURE-SHIFT EXPERIMENTS 

In order to directly probe the influence of the barri- 
ers on the aging dynamics, we perform Jiumerically the 
analog of temperature-shift experimentsEL-,Similar simu- 
lations were already performed in Refs.li-ffi3 with qualita- 
tive results only. Hece we go much further and perform, 
as was done in Ref.til, a detailed quantitative analysis of 
the results. The protocol is the following. The system is 
quenched from T = oo to the temperature T\ — ± AT 
where it ages; AT here will always be positive. At time 
t w , the temperature is shifted to T2, where the measure- 
ments start. We shall discuss the behavior of the auto- 
correlation function C(t w + t, t w ). 



A. Aging is less efficient at lower temperatures 

Let us start with the phenomenology. We observe 
that the decay of the correlation function after the 
temperature-shift from T\ = T2 + AT to T% is slower 
than a purely isothermal aging at T2, meaning that pre- 
liminary aging at a slightly higher temperature 'helps' 
the system at the final temperature. The opposite effect 
is found when T x = T 2 - AT, see Fig. | (top). In that 
sense, aging is less efficient at lower temperatures. 

Moreover, the decay following the shift has the same 
functional form, for small enough AT, as in a simple ag- 
ing experiment at temperature T2. An example of this 
feature is shown in Fig. |[ This implies that the corre- 
lation function after the shift can be superposed to the 
correlation obtained in isothermal aging at T2 by intro- 
ducing an effective waiting time. One has < t w for 
T2 — AT < T2. The same effect is.observed experimen- 
tally when AT is sufficiently smallQ. 

The determination of the effective age of the sample 
can be made rather precise when the results of section [II 
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are used. The correlation function obtained in the shift 
experiment for different AT can be collapsed on the mas- 
ter curves of Fig. ||, using tf^ as a single adjustable pa- 
rameter. This leads to a very precise determination of 
, which does not require any analytical fit of the data, 
see Fig. §, 
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FIG. 8. Comparison between simple aging and shift ex- 
periments in d — 3. Top: Lines are simple aging curves, 
t w = 3728 at T = 0.8 and t w = 1245, 2154, and 3728 at 
T = 0.9 (from bottom to top). Circles are obtained in the 
shift Ti = T 2 - AT = 0.8 -> T 2 = 0.9 at time t w = 3728. One 
concludes that 1245 < tjf < 3728. Bottom: The correlation 
obtained in the shift is superposed to the correlation obtained 



in isothermal aging at T2, which gives the value t% 



2200. 



effective age can be determined by the relation 

<?t 2 -at(U = irAO- (21) 



This relation will be correct if AT is not too large, such 
that quasi-equilibrated structures of sizes < £T 2 -AT(t w ) 
are almost unchanged by the temperature shift. In this 
case, aging at T% is a simple continuation of aging at 
T 2 — AT, at a slightly different rate given by Eq. ([!]). 

Figure || shows that Eq. (12)]), works very well. Such a 
relation was proposed in RefE-3 assuming a pure power 
law growth of the coherence length at all temperatures. 
We show below that our data indeed support Eq. pi]), 
but is incompatible with a pure power law growth of £t- 
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FIG. 9. This curve shows that the effective age of the sys- 
tem is well described by Eq. (^l|). The parameters are the 
same as in Fig. pi 
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B. 'Time is length': Link with the coherence length 



TABLE II. Effective 
shift-experiments . 



waiting times for various 



The effective age of the sample may be simply inter- 
preted in terms of length scales. The growth of the co- 
herence length being slower at lower temperatures, one 
has £T 2 -AT{t w ) < £T 2 {tw)- If one assumes that the age 
of the sample is fully encoded in the value of It, then the 
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Different values of i^f obtained in a series of shift ex- 
periments are reported in Table |l[ This effective age is 
also compared to the simple activation prediction, t^ ct , 
where the same barriers are crossed at the two tempera- 
tures. In this case, one gets the prediction that 



In 



to 



J 2 V T 



(22) 



In this equation, r is the microscopic time that was 
extracted from the growth of the coherence length at 
T = T c . From Table [n| one clearly concludes that 
t~* > i^ ct , which suggests that the microscopic 'trial' 
time is actually much larger than to. The same system- 
atic effect has-been deduced from recent experiments on 
Ising samplesEll 

It is interesting to remark that the simple power law 
growth i T ~ (t/T ) 1/z(T \ with z(T) = z c T c /T, also leads 
to the purely activated law Eq. for the effective wait- 
ing time, and therefore fails to explain the observed be- 
havior of spin glasses. On the other hand, the mixed 
critical/activated growth law described in Eq. ([!]), where 
the microscopic time To is multiplied b y ffi , is indeed 
able to account for deviations-from Eq. (|22|), as already 
discussed in details in Refs.EK3. j-. 

We have used the analysis proposed in Refu to extract 
ip and To from the data in Table |l[ Interestingly, we 
find tjj ~ 1.1, To ~ 2.1, compatible with the value ob- 
tained from the direct fit of the coherence length. The 
agreement between direct and indirect determinations of 
ip is an important result of this paper, since it validates 
the analysis performed on experimental data, where no 
direct determination is possible. However, the value of 
ip = 1.0 favoured by our numerical data is different from 
the ones reported in previous experimental work on Ising 
spin-glasses using different procedures: ip ~ 0.3 — 0.5cj, 
ip ~ 0.7E3, ip ~ It is true that the length scales 

probed in experiments are at least a factor ten larger 
than those probed here. This does not explain, however, 
the scattering of the experimental data. 



VI. LARGE TEMPERATURE SHIFTS: 
REJUVENATION, KOVACS AND MEMORY 
EFFECTS 



A. Is rejuvenation observable in simulations? 

The basic message of large temperature shift experi- 
ments is that, independently of the siga of T\ — T27 aging 
is 'restarted' at the new temperature^. This 'rejuvena- 
tion effect' can be nicely observed through the measure- 
ment of the magnetic susceptibility t w ). For a given 
frequency ui, the dominant contribution to the aging part 
of x(cl>, t w ) comes from the modes with a relaxation time 
~ lo^ 1 which are still out of equilibrium at time t w . Re- 
juvenation after a negative temperature shift comes from 
fast modes, which were equilibrated at Ti, but fall out 
of equilibrium and are slow at Ti- Therefore, one should 
expect to see this phenomenon if the equilibrium confor- 
mation of length scales < It-J^w) is sufficiently different 
at the two temperatures (see below for a more precise 
statement). 

This mechanism is qualitatively different from the pr 
terpretation involving the notion of temperature chaosE-3, 
and put forward in various approachesa&Ell In the 
latter, the existence of an overlap length l ^T\,T%), di- 
verging when T2 — T\ — > 0, is postulated. Its physi- 
cal content is that length scales smaller than l a are es- 
sentially unaffected by a temperature shift Tj — * T 2 , 
while larger length scales are completely re-shuffled by 
the shift. In this picture, rejuvenation is thus attributed 
to large length scales. Strong rejuvenation effects there- 
fore require a very small l Q . 
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We turn how to another set of experimentaJ, where 
larger shiftstl T\ — > T2, and possibly cycles T\ — > T2 — > 
Ti, are performed. In the previous section, indeed, the 
dynamics after a shift was the continuation of the aging 
before the shift. In this section, we use larger tempera- 
ture shifts, so that the small scale structures that were 
equilibrated at the first temperature have to adapt to the 
new one. Precisely how this happens is what we address 
in this section. 



FIG. 10. The 4-point correlation functions at two temper- 
atures in d — 3 (black symbols) and d = 4 (open symbols). 
Times are chosen so that £0.9 ~ ^0.5, namely t(T — 0.9) = 720 
and t(T = 0.5) = 19307. 

It turns out that no clear rejuvenation effects have ever 
been obs erved in simulations of the 3 dimensional Ising 
spin glassli-3E3. This was first attributed to the fact that 
l (Ti, T2) was perhaps numerically large, so that no large 
scale reorganization could-be observed on the time scale 
of numerical simulationaij. Another possibility is that 
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the Edwards-Anderson model-lacks a crucial ingredient 
to reproduce the experiments!!!], or that the length -apd 
time scales involved in the simulations are too smallE3. 

From the above discussion, we see that the crucial in- 
gredient is the small scale reorganization due to a temper- 
ature shift. A natural measure of the spatial organization 
is provided by the 4-point correlation, Eq. (||). We show 
in Fig. |l0|the function C^ir^t) at two different tempera- 
tures Ti = 0.9 and T 2 = 0.5 in d = 3 and d = 4. Times 
are chosen so that It x ~ £t 2 ■ It is clear from this Figure 
that a temperature shift T\ — ► T 2 will hardly play any 
role in d = 3, whereas the two curves are clearly different 
in d = 4. Another way to see this is to observe the tem- 
perature dependence of the exponent a(T) reported in 
Table §. This exponent is almost constant in d = 3, but 
varies significantly in d = 4. This observation suggests 
that no clear effect can be seen in d = 3, whereas d = 4 
should be more favourable. 



B. Negative temperature shifts and rejuvenation 



lo > #Ti(t*)' There should thus exists a typical shift- 
amplitude, AT*(t s ), such that for AT < AT*(t s ), reju- 
vena tion s hould be almost absent. As discussed in sec- 
tion VII B below, we actually can rule out more directly 
this interpretation in terms of an overlap length. 

Although a rather strong rejuvenation appears for 
large AT, one needs to discuss the effect in more de- 
tails. In particular, the experiments show that for large 
AT, rejuvenation is 'complete' in the sense that x( u ,t w ) 
after the temperature shift is indistinguishable from the 
curve obtained after a direct quench from high temper- 
atures. This is closely related to the absence of cooling 
rate effects on the a.c. susceptibility, as reported in Ref.Q. 

We have thus compared the evolution of x(w, t w ) from 
our simulation of a temperature shift to the result ob- 
tained after a direct quench. We find that the curves 
are significantly different. The curve after the shift is 
clearly 'older' than after a direct quench. However, as 
we discuss in the next section, an experimental quench is 
never infinitely fast, contrarily to what can be achieved 
numerically. 



This is indeed what we observe numerically on the ana- 
logue of the a.c. susceptibility, defined in Eq. (||)_Jn 
d = 3, the amplitude of rejuvenation is very smallEj'E3, 
as expected from the behaviour of the 4-point correla- 
tion. In d = 4, on the other hand, the a.c. susceptibility 
'restarts aging' after a negative shift T\ — ► T 2 < Ti, at 
time t s , as illustrated in Fig. 111]. These curves are very 
similar to what is observed experimentally. 
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FIG. 11. Effect of a negative temperature shift on the 
'a.c' correlation function, for ui = 1/1059, after a long stay 
t a = 6449 at Ti = 0.9, for different T 2 , in d = 4. The larger 
AT, the stronger the rejuvenation. 

The amplitude of the rejuvenation is found to increase 
smoothly with the amplitude of the shift AT = T\ — T 2 . 
The 'temperature chaos' picture suggests a more bru- 
tal crossover: no rejuvenation should appear as long as 



C. Cooling rate effects and sub-aging 

In order to quantify the rejuvenation effect, we investi- 
gate the influence of the time t s spent at T\ = 0.9 before 
the temperature shift. The evolution of xi^^w) after 
the shift to T 2 = 0.5 for different t s is shown in Fig. [l^. 
We note that as soon as t s is sufficiently long, t s > 240, 
the evolution after the shift becomes independent of t s . 
For smaller t s , on the other hand, one sees that extra ag- 
ing contributions are present. Hence, for t s > 240, some 
short scale correlations created at T\ survive at T 2 , even 
for large AT, making the relaxation different from what 
it is when t s = 0. This points towards the absence of 
tem peratu re chaos and will be discussed further in sec- 
VII B . This shows also that for large enough t s 



tion 



the 



system behaves after the shift as if it had spent an infinite 
time at T±, i.e. as if t s = oo. 

The important point now is that experiments always 
spend some finite time (actually quite long compared to 
the microscopic time) at all temperatures above the fi- 
nal one T 2 , where some particularly strong correlations 
very rapidly set in and survive when the temperature is 
lowered. Therefore, as soon as the cooling rate is not ex- 
tremely fast, the initial configuration at T 2 already has 
some of the corr elation s that the system wants to grow 
(see also section VIIB). On the other hand, as our sim- 
ulations show, waiting longer at these intermediate tem- 
peratures will not affect further the behaviour at T 2 . The 
initial age of the system is thus effectively non-zero, but 
very soon independent of the cooling rate. 

Interestingly, this non-zero initial age induces apparent 
sub-aging effects. Indeed, if t le \(t w ) = to + t w , where to 
approximately accounts for the aging accumulated on the 
cooling path, the effective exponent n is found to be less 
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than unity: 



dlogt 



rel 



d\og U 



1 + 



< 1. 



(23) 



We have confirmed this directly on the scaling of the two- 
time correlation function C(t + t w ,t w ) obtained after a 
slow quench in the d = 4 case. We find that /i = 0.96 < 
1, whereas the scaling obtained after an infinitely fast 
quench indicated super-aging, jj, = 1.05 > 1 (see Fig. |] 
above). We believe that this effect is significant. It is 
thus tempting to ascribe at least part of the sub-aging 
effects seen experimentally to finite cooling rate effects. 
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FIG. 12. Effect of a negative temperature shift on the 'a.c' 
correlation function, after a stay at Ti = 0.9 of various dura- 
tions t s . Final temperature is T2 = 0.4. 



length scales as compared to the length scales involved in 
the aging at T\. Thus, when the temperature is shifted 
back to Ti, the correlations of length scale £t 2 (t' s ) grown 
at T2 almost instantaneously re-equilibrate at T± (in fact 
in a 'memory' time scale t m such that ^Ti(im) ~ ^T 2 {t' s )^ 
which implies that t m <C t' a when AT is sufficiently 
large). The memory is thus stored in the intermediate 
length scales, between #r 2 (0 an d ^Ti(i s )- 

Hence, the explanation of the memory effect relies on 
the separation of length scales only. This ingredient is 
distinct from the one needed to observe rejuvenation, 
which relics on the reorganization of small length scales 
after a temperature change. Since we have shown that 
these two ingredients are present in the 4 dimensional 
spin glass, we are able to reproduce experimental data 
very well in Fig. where, for purely esthetic reasons, a 
double cycle was performed. 
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D. Temperature cycles and memory 

Since we have considered above both cases of a posi- 
tive and a negative temperature cycle, we are now in posi- 
tion to combine both procedures and study temperature- 
cycles. The experimental procedure is here T = 00 — > 
T\ — ► T% < Ti — ► Ti. The time spend at Ti is t s and 
the time spent at T2 is t' s . The spectacular 'memory ef- 
fect' arises when the temperature is shifted back to Ti. 
It is observed that although aging was fully restarted at 
T2, the system has a strong memory of the previous ag- 
ing at Ti . The dynamics at Ti rproceeds almost as if no 
cycle to T2 had been performedu. The coexistence of re- 
juvenation and memory was made more spectacular in 
the 'dip-experiment' proposed in Rei.0. This protocol is 
too complicated to be studied theoretically, but basically 
carries the same physical content as the cycle we discuss 

As discussed i n!0 the memory effect is a simple conse- 
quence of the separation of time and length scales. When 
the system is at T2 < Ti, rejuvenation involves very small 



FIG. 13. Evolution of the 'a.c' correlation function in the 
procedure T = 00 -> T ± = 0.9 -> T 2 = 0.4 -» T ± -» T 2 , 
showing, as in experiments the coexistence of rejuvenation 
and memory effects. 



E. Positive temperature shifts 

From the results of the previous sections, the physics in 
a shift experiment is the following. At the first tempera- 
ture Ti , the system evolves towards equilibrium through 
the growth of a coherence length It x . When the tem- 
perature is shifted to T2 at time t s , all length scales are 
driven out of equilibrium. Length scales smaller than 
^Ti(is) undergo a 'quench' from Ti to T2, while larger 
length scales which were not equilibrated at Ti undergo 
a quench from T = 00 to T 2 . If T 2 < Ti, then larger 
length scales do not matter due to the huge separation 
of time scales. 

The situation is different in a shift such that T% > T±. 
Then small length scales have to 'unfreeze' to find their 
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new equilibrium at T2, while larger ones which where 
frozen at Ti grow as if the quench had been from T = 00. 

To support further this physical picture, we performed 
a positive temperature shift experiment T± = 0.5 — > T% = 
0.9, after time t s — 19307 at Ti, and then some extra time 
t at T2. The results are described in Fig. ^ which shows 
both the behavior of the autocorrelation and the 4-point 
correlation after the shift. 



t = 2 
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6449 o 
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Then, small scale structures very rapidly equilibrate 
at the new temperature, t < 416. Note that it took 
t s = 19307 to reach the same coherence length at Ti, 
a consequence of the length scale separation. After this 
short transient, dynamics proceeds as if the initial stay at 
Ti was not present, and the subsequent aging is very sim- 
ilar to isothermal aging at T2, as soon as t > 2154. The 
same features are also clearly visible on the 4-point corre- 
lation function, see Fig. [l4] (bottom). Note in particular 
how small scales rapidly 'unfreeze' before the large scales 
evolve towards equilibrium: the correlation for t — 240 
is below the one for t = 0, before the coherence length 
£r 3 (i) grows beyond £ Tl (t s ). 



F. The 'Kovacs effect' 

This dual behavior between small and large length 
scales results in a spectacular effect, which was first ob- 
served by Kovacs in polymeric glassestia. Further devel- 
opments may be found in Refs.EjO. Since it is referred 
to in the literature as a 'memory effect', but is different 
from what the spin glass literature names 'memory' (see 
above), we shall follow RefO and describe this as the 
'Kovacs effect'. 




r 



FIG. 14. Autocorrelation function (top) and spa- 
tial correlation function (bottom) in a positive cycle 
Ti = 0.5 — > T2 = 0.9. Times t refer to time spent after 
the shift. For comparison, full lines show curves obtained in 
a direct quench to T2. Small and large scales behave very dif- 
ferently. Note also the similarity between both figures which 
makes clear the link between time and length. 

Let us comment first on the time correlation functions, 
Fig. [l4| (top). Immediately after the shift, t = 0, the 
decay of the correlation shows a short-time part which 
is slower than the reference curve with t w = 19307 at 
T2, and a long-time part which is faster. This nicely 
illustrates the two types of structures present at that time 
in the system. 
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FIG. 15. Behavior of the energy density in a shift experi- 
ment from Tl = 0.6 to T2 = 0.7, at time t w — 3728, compared 
with its decay in simple aging experiments (lines), in d = 3. 
A simple shift with t^f = 1450 allows one to superimpose the 
curves. 

We focus here on the energy density e(t) after a pos- 
itive temperature shift. The Kovacs effect concerns the 
specific volume, but the difference is irrelevant for our 
purposes. Like in recent numerical experimentaL-3, we 
find that the decay of the energy density following the 
shift follows the same time evolution as in the simple ag- 
ing case, if an appropriate effective waiting time ^ is 
properly taken into account, see Fig. |lj. We find that 
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the effective age of the sample defined from the correla- 
tion function or from Eq. ( |2l| ) works well for the energy 
density also. This is illustrated in Fig. [l5| . 

KovacsE3 noticed that in a similar protocol on polymer 
glasses, the same non-monotonic initial behavior could be 
seen in the evolution of the specific volume as we observe 
in Fig. [ITJ for the energy density of the spin glass imme- 
diately after the shift. Zooming on the transient region 
and setting the origin of time when the temperature is 
shifted leads to the curves plotted in Fig. |l6|. The top 
curves of Fig. 16 are specially designed to follow Kovacs' 
experiments, where the time t s of the shift is chosen so 
that eTi(is) = e.T 2 {t = 00). Since the energy density has 
already the correct equilibrium value at the new temper- 
ature, the naive expectation is that e(t > t s ) = const. 
Instead, the non-monotonic behavior of Figs. [l6| is ob- 
served. 
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FIG. 16. Non-monotonic behavior of the energy in a shift 
experiment. Top: Same T2 but different AT. Waiting times 
are such that the energy just before the shift is ~ e cq (T2) as 
in experiments on polymeric glasses. Bottom: The AT is the 
same in all shifts, but shift times are different, t 3 = 57797, 
19307, 6449, 2154, 720 and t s = (from bottom to top). 



The presence of a growing coherence length allows one 
to gitte a very simple interpretation of this 'Kovacs ef- 
fect 'ta. It results indeed precisely from the dual behavior 
of length scales described above. When the temperature 
is shifted to T2, length scales shorter than trAts) have 
to re-equilibrate at T2, where their equilibrium energy is 
higher than at T\. This explains the initial rise of e(t s +t). 
On the other hand, length scales larger than i^its) still 
have to 'cool down' and decrease their energy. These two 
opposite trends directly explain the 'Kovacs effect'. This 
scenario was recently illustrated on the exactly soluble 
example the 2D XY modeM 

It is possible to be more quantitative here, using 
the coherence length as an ingredientta. The time 
scale tx where the energy density reaches its maxi- 
mum corresponds in this picture to the time where small 
length scales have re-equilibrated at the new tempera- 
ture. Hence, an excellent approximation for this time 
scale should be: 



£ T2 (t K ) ~ l Tl (t s ). 



(24) 



This relation says that tx is an increasing function of t s , 
and a decreasing function of the temperature difference, 
as is obvious from Figs. [16[ The height ex of Kovacs' 
hump varies in the opposite direction, as expected from 
the inverse power-law dependence of the excess energy 
with— the coherence length, ex ~ £T 2 (tK) e ~ d , found in 
Ref.EJ. We numerically find that all the curves of Figs. [l6| 
can actually be collapsed into a single master-curve (see 
Fig. [I?]) which thus takes the form: 



e(t + t s ) ~ e oq (T 2 ) + e K f 



(£)• 



(25) 



where the scaling function f(x) behaves approximately 
as a power law for large arguments. 
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FIG. 17. All the curves of Fig. [j~6| are superposed using 
Eq. (Eq). The curve 't s = 0' refers to a direct quench to T2. 
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VII. PHYSICAL DISCUSSION 

A. Physical picture of the spin-glass phase 

Mean-field theories have nothing to say about possi- 
ble relevant length scales (such as £t) and their time and 
temperature dependence. On the other hand, the droplet 
picture, which focuses on relevant length scales, seems to 
miss some important points such as the power-law be- 
haviour of the 4-point correlation function. This point is 
important since it allowed us to account for rejuvenation 
effects without the need of the concept of temperature 
chaos. 

One interpretation of the above results is that, as pre- 
dicted by mean-field theories and given some credi t by re- 
cent numerical work on low-lying excitations" 6 H 17 1, differ- 
ent equilibrium configurations are accessible to the spin- 
glass in its low temperature phase. These configurations 
have a global overlap which is close to zero, but can be lo- 
cally similar. The fact that the stationary part of C^r, t) 
decays as a power-law suggest the existence of a fractal 
'backbone' of spins that have identical mutual orienta- 
tions for all these configurations, with a fractal dimen- 
sion df — d— a. It is reasonable to assume that the small 
scale properties of this backbone will be temperature de- 
pendent: more spins will freeze and join the backbone as 
the temperature is reduced. The simplest scenario com- 
patible with a zero minimal overlap is that the backbone 
is dense on small scales, and fractal on large scales, with 
a temperature dependent crossover length I* (T) . The ef- 
fective exponent a(T) would in this case decrease with 
temperature, as seen numerically. Another possibility is 
that the fractal dimension (and thus the exponent a) is 
truly temperature dependent, as in the low temperature 
phase of the 2D XY model. As discussed above, there are 
actually many phenomenological similarities between the 
spin-glass phase and the 2D XY model, provided length 
scales, rather than time scales, are comparedEjeH. 

The power-law decay of C±{r, t) suggests that the 
whole spin-glass phase is in a certain sense critical, at 
least in the 'zero-overlap' sector which aaas indeed found 
to be mass-less in replica field analysisu. However, this 
is not in contradiction with the existence of a finite corre- 
lation length £(T) separating critical from activated dy- 
namics within a single 'state'. A pure power-law growth 
of It is not necessarily a consequence of the criticality of 
the spin-glass phase. 

Of course, the numerical evidence for this scenario is 
fragile, and it could be that C±(r,t — > oo) in fact tends 
for large r towards gj| A . For the purpose of interpreting 
aging experiments, however, it is sufficient that this sce- 
nario holds even approximately on the relevant time and 
length scales. 



B. Rejuvenation from small scales and absence of 
temperature chaos 

The role of temperature changes can be exactly com- 
puted in the Random Energy ModelEil and in the critical 
phase of the XY modeO. Both show that it is possi- 
ble to induce strong rejuvenation effect without the ex- 
istence of an overlap length. Several facts, reviewed in 
Ref.M, also suggest that the overlap length is not rel- 
evant to the experimental findings. Here, we want to 
address this question more precisely on the basis of nu- 
merical results. As shown above, we can now observe 
beyond any doubts rejuvenation (and memory) effects 
in the d — 4 Edwards- Anderson model which are very 
similar to those observed experimentally. We have also 
investigated directly the way configurations evolve dur- 
ing a temperature shift using a mixed 4-point correlation 
function, defined as follows: 



1 



N 

C i (r,£,T 1 ,AT) = -Y,' 

i=l 



s a S w )s a l+r {t w )s\(V w ) S \ +r (t' w )), 



(26) 



where replica a is at temperature T\ , replica b at temper- 
ature T-2 — Ti — AT, and the times t w , t' w are chosen such 
that the coherence length is equal to a common value i 
at the two temperatures. Obviously, when AT = 0, this 
correlation function is identical to the previous one. For 
AT > 0, this correlation function measures the similarity 
between the patterns grown at the two different temper- 
atures. In a temperature chaos scenario, one expects the 
following inequality: 



C 4 (r, i, Ti , AT) < d (r, i,T\, AT = 0) . 



(27) 



10 u ^ r 







10" 



10" 



10" 



- 1 1 1 

T = 0.5, = 19307 -a- 
T= 0.9,i„ = 720 - -a- 
AT = 0.4 -»■ 
AJ = 1.0 
AJ = 0.5 
AJ = 0.1 — r- 



1 



FIG. 18. Mixed correlation functions, both for temperature 
changes and for coupling changes (black symbols), compared 
to standard four-point functions (open symbols). It is clear 
that the two perturbations (AT, AJ) have qualitatively dif- 
ferent effects. 
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Figure [18] shows that this is not the case. The results 
are actually compatible with the idea that the same pat- 
terns grow at the two temperatures - the backbone sup- 
porting the common parts of these patterns being more 
fluffy at lower temperatures. This conclusion was already 
reached above when we discussed cooling rate effects. 

An interesting comparison can be made with a situ- 
ation where-ehaos is expected, e.g. when the couplings 
are changedo. We therefore also show in the same figure 
the mixed correlation function when the couplings are 
changed between replica a and replica b according to: 



J i 



VI + A J 2 



(28) 



where AJjj are independent Gaussian variable of vari- 
ance A J 2 and mean 0. In this case, the inequality (^7|) is 
indeed clearly observed. We conclude thus that AT and 
A J have a qualitatively different influence on the system. 

Note that our mixed correlation function, once inte- 
grated over space, leads to the overlap between the two 
tempecatures. The latter quantity was studied directly 
in Ref.E3 in d = 3, with conclusions similar to ours. 

The simultaneous observation of rejuvenation and ab- 
sence of temperature chaos is an important result of this 
paper. In d = 3, no temperature chaos was found, but no 
rejuvenation either. This left the door open to the pos- 
sibility that the length scales investigated were too small 
to observe these two effects. We have thus demonstrated 
that both issues can be separated. Of course tempera- 
ture chaos on large length scales is still possible, but is 
not needed to interpret rejuvenation effects. 

In summary, our results confirm that rejuvenation is 
due to the freezing of small length scale modes which were 
'fast' at the higher temperature. This freezing changes 
the correlations on small scales, as seen on the 4-point 
correlation function. This is in agreement with the sce- 
nariaJiased on a hierarchy of length scales proposedJxi 
Rcf .BEJ, and with the phenomenology of the XY modeO, 
and is markedly different from the temperature chaos 
picturp- This feature can be illustrated in the 2D XY 
modelEJ, where each Fourier mode ip(q) of the order pa- 
rameter is affected by a temperature shift. In the spin- 
wave approximation, one has (<p(q)) ~ T/q 2 . Hence, each 
mode is affected when the temperature is changed by AT 
by an amount 



AT 

5{<p(q)) ~ —, 

q 2 



(29) 



which shows that larger length scales are more influenced, 
but with no typical 'overlap length'. 



VIII. SUMMARY AND CONCLUSION 

The interest of a long paper is that a detailed discus- 
sion of rather subtle points can be given. The drawback, 
obviously, is that the message is somewhat diluted. We 



therefore give in this last section the main conclusions 
from our study and end on open problems. 

• Aging dynamics in spin-glasses can be associated 
with the growth of a coherence length £t, sepa- 
rating small, equilibrated scales < £t from large 
frozen, out of equilibrium scales > £t- This scale 
is however not a domain size in the usual sense, 
but rather the size of a backbone of spins common 
to all spin-glass configurations. This interpretation 
stems from the power-law decay of the 4-point cor- 
relation function from which £t is extracted. 

• The coherence length £t follows a critical power- 
law growth at small times that becomes activated 
for larger times, and is well described by Eq. (ffl). 
The associated barriers T(T) vanish at the criti- 
cal temperature. The barrier exponent ip was es- 
timated to be -0 ~ 1-0 for d = 3 and ip ~ 2.3 in 
d = A. 

• This mixed critical/activated growth law allows one 
to interpret several important aspects of both simu- 
lations and experiments, for example the deviations 
from a purely activated behaviour that are revealed 
by temperature shift procedures, or the super-aging 
behaviour of the correlation function observed in 
d = 4. 

• The short scale behaviour of the 4-point correla- 
tion is quite sensitive to temperature in d = 4, but 
much less in d = 3. This in turn leads to strong 
rejuvenation effects in d — 4, quite similar to those 
observed in experiments, that we observe for the 
first time in simulations. 

• An interpretation of the observed rejuvenation in 
terms of temperature chaos is, we believe, ruled 
out: see Fig. JIs|). Rather, some correlations built 
at a higher temperature persist and are reinforced 
at lower temperatures. 

• A finite cooling rate effect follows from this, which, 
interestingly, leads to an apparent sub-aging be- 
haviour for the correlation function, instead of the 
super-aging that holds for an infinitely fast quench. 
The cooling rate dependence however saturates 
quickly as soon as the cooling rate is not infinitely 
fast. Both these features agree with experiments, 
for which the cooling rate is always very slow com- 
pared to microscopic frequencies. 

• The dichotomy between small, equilibrated scales 
and large, frozen scales allows one to account semi- 
quantitatively for many features, such as the role 
of temperature shifts, the memory effect or the Ko- 
vacs' hump. 

Although our results are suggestive, several unsettled 
points remain. In particular, rejuvenation effects are 
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found in d = 4, but not in d = 3, whereas experiments are 
obviously performed in d = 3. We conjecture that for the 
time scales investigated, the large scale topology of space 
is irrelevant, and the major difference between d = 3 and 
d = 4 should rather come from the local connectivity. 
Hence it should be possible to obtain rejuvenation in a 
d = 3 model with more neighbours, and reproduce most 
experimental results with a realistic model. 

The most important theoretical point is obviously the 
nature of the spin-glass phase. A well posed problem 
(but very difficult to settle numerically) is the true long 
distance behaviour of the 4-point correlation function: 
power-law decay, as expected from replica symmetry 
breaking theories, or convergence towards q% A , as for a 
disguised ferromagnet? The final picture of real spin- 
glasses might in the end have borrow concepts from both 
theories. The hope is that the concepts that will emerge 
will be useful to understand many other glassy systems, 
which share a very similar phenomenology. 
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NOTE ADDED 

After this manuscript appeared as a preprint (cond 
mat /0202069 ), a paper by Yoshino et al. 



4 E. Vincent, J. Hammann, M. Ocio, J.-P. Bouchaud and 
L. F. Cugliandolo, in Complex behavior of glassy systems, 
Ed.: M. Rubi Springer Verlag Lecture Notes in P hysics, 



( cond- 



mat / 0203267| ) appeared where the dynamics of 4-d EA 
model is studied, and the results also interpreted as a 
crossover between critical and activated dynamics. The 
value of the exponents z c and ip given in that paper 
slightly differ from those obtained here. For example, 
ip is found to be in the range 2.5 — 3 whereas we report 
■0 ~ 2.3. One possible explanation is that the procedure 
to extract ^r(ttu) form Ci(r,t w ) is quite different. 
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